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Abstract 

A two dimensional self-gravitating Hamiltonian model made by N fully- 
coupled classical particles exhibits a transition from a collapsing phase (CP) 
at low energy to a homogeneous phase (HP) at high energy. From a dynamical 
point of view, the two phases are characterized by two distinct single-particle 
motions : namely, super diffusive in the CP and ballistic in the HP. Anomalous 
diffusion is observed up to a time r that increases linearly with N. There- 
fore, the finite particle number acts like a white noise source for the system, 

inhibiting anomalous transport at longer times. 
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In the past years the thermodynamical properties of gravitatial models have been studied 
in detail from a theoretical and computational |2|-|J point of view. In particular, it has 
been shown that at low energy the gravitational forces give rise to a collapsing phase (CP), 
identified by the presence of a single cluster of particles floating in a diluted homogeneous 
background. At high energy a homogeneous phase (HP) is recovered: the cluster disappears 
and the particles move almost freely. In the transition region the system is characterized 
(in the microcanonical ensemble) by a negative specific heat: the corresponding instability 
(termed "gravo-thermal catastrophe") is of extreme relevance for astrophysics (see || for 
more details). This apparent thermodynamical inconsistency has been solved by Hertel 
and Thirring in Ref. 0, where they demonstrated the non equivalence of canonical and 
microcanonical ensemble in proximity of the transition region. These theoretical results 
have been successfully confirmed by numerical investigations of self-gravitating non-singular 
systems with short range interaction 0|§. 

More recently, in one dimensional lattices of fully and nearest-neighbour coupled sym- 
plectic maps with an attractive interaction it has been noticed that clustering phenomena 
are associated with anomalous diffusion (in particular, with subdiffusive motion), at least 
for short times ||. Anomalous diffusion can be defined through the time dependence of the 
single particle mean square displacement (MSQD) < r 2 (t) >, that typically reads as 

<r 2 (t)>ott a (1) 

where the average < • > is performed over different time origins and over all the particles 
of the system. The transport is anomalous when a ^ 1: superdiffusive for 1 < a < 2 and 
subdiffusive for < a < 1 |7]||. Anomalous transport has been revealed in dissipative and 
Hamiltonian models J7| as well as in experimental measurements ||. However, the main 
part of literature focuses on systems with few degrees of freedom (namely, one or two) and 
only few studies have been devoted to extended models with N ^> 1 [[]]. 

In this Rapid Communication, the thermodynamical and dynamical properties of a 2- 
D Hamiltonian system, constituted of N particles interacting via a long range attractive 
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potential, are analyzed. In particular, we observe a transition from CP to HP associated 
to a dynamical transition from anomalous to ballistic transport. Finite N effects induce 
a crossover from anomalous to normal diffusion at long times. In the limit N — > oo, the 
transport mechanism remains anomalous for any time and reduces to that of a single particle 
in an "egg-crate" potential |TD| . 



We consider a system of N identical fully coupled particles with unitary mass evolving 
in a 2 dimensional periodic cell described by the Hamiltonian fllllJl2 



n 2 + 2 , N 

H = K + V = y Px >* Py ' 1 + —y 

^ 2 2N 4i 



cos(xi - Xj) - cos(yi -yj) - cos(xi - xj) cos(t/, - yj) 

(2) 



where (xi,p Xji ) and (yi,p yt i) are the two pairs of conjugate variables with {x^yi) G 
[— 7r, 7r[x [— 7i, 7r[, K and V are the kinetic and potential energy, respectively. The poten- 
tial part corresponds to the first three terms of the Fourier expansion of a 2-D attractive 
potential of the kind V(r) oc log | r | . Such type of interaction arises in self-gravitating 2- 
d gases [|l]-[J as well as in point vortices model for 2-d turbulence [13]. Due to the long 



range interaction among all the particles, this model can be described in terms of meanfield 
variables. In particular, the potential energy can be rewritten as V = | Y^iLi V%i with 



Vi = 3- M x cos(xi - <j> x ) - My cos(yi - (j) y ) - ]- \M xy cos(xi + y { - <j>+ y ) + M xy cos(xi - y { rll . 



(3) 

where M z = (< cos(z) >at,< sin(2;) > N ) = M z exp[i(j) z ] represents four two-dimensional 
meanfield vectors with z = x,y,x ± y and < .. >n denotes the average over N. However, 
the single particle potentials Vi are non-autonomous since the meanfield quantities M z and 
<p z are defined through the instantaneous values of the particles coordinates. The motion of 
each particle is therefore determined self-consistently by an attractive and non-autonomous 
force field that is itself uniquely determined throught the motion of all the particles. The 
self-gravitating nature of the model is due to this effective force acting among the particles. 
Since V is invariant under the transformations x ^ —x, y <-> —y and x <-> y, it turns out 



that in the "meanfield limit" (i.e. for N -> oo with U = H/N constant) M x = M y = M 
and M+y = M~ y = P. Moreover, in this limit and assuming <fr z = the single potential V% 
turns out to be an egg-crate potential similar to that studied in JTIJ . This periodic potential 
is characterized in each elementary cell by a minimum (V m = 3 — 2M — P), 4 maxima 
(V M = 3 + 2M - P) and 4 saddle points (V a = 3 + P). 

For specific energy U smaller than a critical value U c ~ 2, we observe that the particles 
are mainly in a clustered state. Above U c a HP is recovered. Following Refs. |6j,fL^|, the 
degree of clustering of the particles can be characterized through the time averages < M z > t 



14|,[15| . When at each time the particles have almost the same position < M z > t are 0(1), 



while for a HP their values vanishes as l/y/N |12| . Fig. 1 shows that, for U — ► 0, the average 



quantities M, P tend to one. This indicates that the particles are almost all trapped in a 
potential well of depth ~ (V s — V m ) forming a compact cluster. For increasing energy U, the 
kinetic contribution becomes more relevant and the average number of particles trapped in 
the potential well reduces. As a consequence the value of < M z > t decreases together with 
(V 8 — V m ). For U > U c , the system is no more clustered and the particles can move almost 
freely. Moreover, due to finite N effects < M z > t is not exactly zero, but 0(1/ yN). 

In Fig. 2 the temperature T =< K > t /N is reported as a function of U. Above U c , T 
increases linearly with U indicating that the system behaves like a free particle gas. In the 
CP, the tendency of the system to collapse is balanced by the increase of the kinetic energy 
0. This competition leads initially (for < U < 1.8) to a steady increase of T, followed (for 
1.8 < U < U c ) by a rapid decay of T. This yields a negative specific heat as illustrated in the 
inset of Fig. 2. These results are in full agreement with theoretical predictions based on the 
analysis of a simple classical cell model and with numerical findings for short ranged 
attractive potentials. The phenomenon of negative specific heat can be explained within 
a microcanonical approach with an heuristic argument |1[: approaching the transition, a 
small increase of U leads to a significative reduction of the number of collapsed particles (as 
confirmed from the drop exhibited by M and P for U > 1.8); as a consequence the value of 
V grows and, due to energy conservation, the system becomes cooler. 



Our data confirm also another important prediction of Hertel and Thirring |lj: the non- 
equivalence of canonical and microcanonical ensemble nearby the transition region. In the 
inset of Fig. 2 are reported the microcanonical findings, obtained via standard molecular 
dynamics (MD) simulations, and the theoretical canonical results, derived in the mean-field 
limit These two sets of data coincide everywhere, except in the energy interval 

1.6 < U < 2.0. The discrepancy is due to the impossibility of the canonical ensemble to 
exhibit a negative specific heat, a prohibition that does not hold for the microcanonical 
ensemble. Our theoretical estimation of the Helmholtz free energy F = F(T) reveals that 
usually F has an unique minimum. For T < 0.5 the minimum Fq corresponds to non zero 
values of M and P (i.e. to the CP), while for T > 0.55 the minimum Fh is associated to 
M = P = (i.e. to the HP). In the region 0.5 < T < 0.55, both minima Fc and Fh coexist 
as local minima of the free energy. However, for T < T c = 0.54 the CP is observed because 
Fc < Fh , while for T > T c the HP prevails since Fh < Fc- At T = T c the two minima 
are equivalent and a jump in energy from U(T~) ~ 1.6 to £/(T c + ) « 2.0 is observed. This 
picture suggests that this transition can be considered as a first order transition 0. 

Let us now investigate if the observed thermodynamical transition has any consequence 
on the dynamical behaviour of the system. In order to characterize the single particle 
dynamics, we consider the MSQD < r 2 (t) >. As shown in Fig. 3, in the CP the diffusion is 
anomalous for times shorter than a crossover time r, while for longer times the Einstein law 
is recovered < r 2 (t) >oc 4Dt (where D is the diffusion coefficient). A similar behaviour for 
the MSQD has been already observed for a system of N coupled symplectic maps in 0, but 
with a < 1. However, in the present case r increases linearly with N |15[ indicating that in 



the meanfield limit the asymptotic dynamical regime will be superdiffusive JTH 



The observed dynamical behaviour can be explained noticing that in the mean-field limit 
each particle i will see essentially the same constant 2-D egg-crate potential V{. Moreover, 
it has been shown in Ref. |T0| that a single particle moving in a egg-crate potential with an 



energy between V s and Vm exhibits super diffusion. This is due to the fact that the particle 
moves for long times almost freely along the channels of the potential and episodically is 



trapped for a while in the potential well. In phase space, the super diffusive phenomenon 
can be explained by a trapping mechanism in a hierarchy of cantori around a cylindrical 
KAM-surface 0. Therefore, in our model (|2|) for N —>■ oo anomalous transport is due only 
to the fraction of particles that can move along the channels. 

For finite N, the potential Vi seen by the particle % will fluctuate in time. Hence, particles 
having an energy close to Vs have the possibility to be trapped in the potential well as well 
as to escape from it. As a consequence, for sufficiently long time scales each particle can 
experience free and localized motions. The fluctuations of the potential Vi reflect themselves 
on the structure of the phase space, introducing a white noise that destroys the self-similar 
structure of the island chains and of the cantori below a certain cut-off size. Being the 
self-similarity no more complete, one expect that on long time scales normal diffusion will 
be recovered fll9| . 

As pointed out in Refs. pOplf , if white noise is added to a dynamical system exhibiting 
super diffusive behaviour, D (measured in the limit t — > oo) is inversely proportional to the 
noise amplitude. Therefore, we expect that in our model the value of D will increase with 
N. That is indeed the case, and we observe a power-law dependence of the type D oc iV 7 . 
For example, considering systems with 100 < N < 10,000 we have found for U = 1.48 and 
U = 2.00 a 7- value equal to 0.7 ± 0.1 and 1.0 ± 0.1, respectively. The iV-dependence of the 
diffusion coefficient can be explained noticing that D oc r a_1 |f?2|| . This result coincide with 
that found theoretically in Ref. |2(J and confirmed numerically by considering as dynamical 



models two very simple noisy maps. For subdiffusive motion D is inversely proportional to r 
(as found in [[HJ), while for super diffusive motion (a > 1) a direct proportionality is expected 
p0| . As already reported r oc N, therefore we will have that 7 = a — 1. Assuming for a the 



corresponding asymptotic values P3| , we can estimate as theoretical values 7 ~ 0.64 and 
~ 0.9 for U = 1.48 and 2.00, respectively. In view of all the present limitations, these values 
can be considered consistent with the numerical findings. 

As a final point, we examine the dependence of the asymptotic a- values from the energy 
U of the system. A transition from anomalous diffusion to ballistic motion (a = 2) at 
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U ~ U c is evident from Fig. 1, where the a- values, obtained for N = 4, 000, are reported. 
In particular, for 0.4 < U < 2.0, we observe an increase of a from 1.3 ± 0.1 to 1.9 ± 0.1. 
This phenomenon is a consequence of the flattening of the single particle potential (i.e. of 
the reduction of Vm — V m ) observed for growing U. The decrease of the average number of 
particles trapped in the cluster, and the consequent increase of those moving freely, naturally 
drives the diffusion mechanism toward a ballistic behaviour. Moreover, for U > U c , the 
potential Vi fluctuates with typical amplitude 0(1/ y/N) around a constant value and a 
ballistic motion is expected for all the particles. For small energies (U < 0.3) the MSQD 
seems to saturate to a constant value, indicating that all the particles are always clustered. 
However, we believe that in these cases our observation time was not sufficient to detect 
particles escaping from the potential well. 

In conclusion, we have shown for the first time a thermodynamical transition associated 
to a dynamical transition from anomalous to ballistic transport. Moreover, the transport 
in our iV-body system can be interpreted in terms of a noisy single particle motion in a 
2-D Hamiltonian egg-crate potential. The asymptotic dynamics of the model is strongly 
influenced by the order in which the two limits N —* oo and t — > oo are taken. Indeed, if the 
limit iV — > oo is performed before the limit t — > oo the diffusion will be always anomalous. 
Otherwise normal diffusion is recovered for sufficiently long times. 

As noticed in Q, the dimensionality of the system should not affect the main character- 
istics of the observed thermodynamical transition. We expect that the same should be true 
also for the corresponding dynamical behaviours. Moreover, anomalous diffusion should be 
observable for atomic clusters |24[] , turbulent vortices |13| (for which indeed has been already 



observed ||) and gravitational systems 0. All the cited systems share as a common aspect 
to exhibit a clustered phase. 

As a final remark, we claim that the inclusion of other terms of the Fourier expansion in 
the expression of the gravitational potential, will not affect the main results here presented. 
However, as far as transport properties are concerned, we believe the range of the force to 
play an important role. In particular, for short-ranged interactions we expect a more chaotic 



behaviour of the particles. This may prevent the anomalous diffusion regime to persist in the 
thermodynamical limit ]nj . Future work will be devoted to the study of these fundamental 
points |T5f] . 
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FIGURES 



FIG. 1. Time averages of M and P as a function of U. The solid curves refer to the theoretical 
estimation (i.e. to canonical results) and the symbols to the MD findings (i.e. to microcanonical 
results). The exponents a, defined in eq. (Q), are also reported (triangles). The MD data have 
been obtained with N = 4, 000 (apart few points with N = 10, 000) and averaged over a total 
integration time ranging from t = 1,2 x 10 6 to t = 2,4 x 10 6 , with a time step dt = 0.3. The 
a-values have been estimated in the time interval 150 < t < 10, 000 for any reported U. 

FIG. 2. Temperature T as a function of the specific energy U. The solid line corresponds to 
the analytical estimation (canonical ensemble) and the triangles to the simulations results (micro- 
canonical ensemble). In the inset, an enlargment of the transition region is reported: the solid 
(resp. dashed) curves refer to the principal (resp. relative) minimum of F(T). The parameters for 
the MD simulations are the same as in Fig. 1. 

FIG. 3. Mean square displacement < r 2 (i) > as a function of time in a log-log plot. The 
cross-over time r from anomalous to normal diffusion is also reported. The data refer to U = 1.1 
and N = 4,000, the total integration time is t = 1.2 x 10 6 . In the inset, the logarithm of the 
cross-over time r is reported as a function of ln(iV) for U = 1.48. The reported r-values (circles) 
have been estimated adopting a threshold (3 = 1.1 (for more details see [17]). The solid line 
represents a best linear fit to the data and its slope is 0.95 ± 0.08. 
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